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ABSTRACT 

We present a systematic analysis of the X-ray spectra of NGC 1313 X-l and NGC 1313 X-2, 
using three years of XMM-Newton observations. We fitted the continuum with a Comptoniza- 
tion model plus a multicolor blackbody disc, that describes the effects of an accretion disc 
plus a corona. We checked the consistency of this spectral model on the basis of the variabil- 
ity patterns of its spectral parameters. We found that the two sources show different spectral 
states. We tentatively interpret the observed behaviour of NGC 1313 X-l and X-2 within the 
framework of near Eddington and/or super-Eddington accretion. We also attempted to deter- 
mine the chemical abundances in the local environment of NGC 1313 X-l and X-2 from the 
EPIC and RGS spectra. The results appear to indicate subsolar metallicity for both sources. 
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1 INTRODUCTION 

Ultra luminous X-ray sources (ULXs) are extragalactic, point-like, 
off nuclear X-ray sources with observed isotropic X-ray luminosi- 
ties higher than the Eddington luminosity for accretion onto a 10 
Mq black hole, and typic ally in the rang e from 10 39 erg s _1 up 
to 



10 4i erg s" 1 (e.g. 



Fabbiano 



1989). Several pieces of ob- 



servational evidence strongly suggest that the majority of ULXs 
are accreting black hole (BH) X-ray binaries (XRBs) with massive 
donors (e.g. Zampieri & Roberts 2009 and references therein). 

At present, the size and the origin of the BHs hosted in ULXs 
is still a matter of debate. It has been proposed that ULXs are 
powered by a stellar mass BHs (10-^20 Mq) accreting largely 
above the Eddington limit and/or having beamed emission (e.g. 



King et al. 200 1[ [King 


2009; Begelman 2006, Socrates & Davis 


2006| Poutanen et al. 


2007). Alternatively, the compact object 



could simply be bigger, and the accretion would be in the usual 
sub-Eddington regime. In this case the compact object would be 
an intermediate mass black hole (IMBH) with a mass in excess 
of 100M Q (e.g. |Colbert & Mushotzky|1999| >. Another interpreta- 
tion follows from the possibility of forming massive stellar remnant 
BHs (20-^80 A/0) from the direct collapse of low-metallicity, mas- 
sive (> 30 - 40 M ) stars ( |Fryer|1999HBelczynski et al.|2010| >. 
Some ULXs may contain such massive stellar remnant BHs accret- 
ing at or slightly above the Eddington rate, formed in low metal- 
licity environments (Z ~ 0. 1 / . : e.g Zampieri & Roberts 2009 
|Mapelli et al.|2009||2"0"T0l >. 

The spectra of many ULXs are qualitatively similar to those of 
Galactic XRBs, and can be well reproduced by a multicolour disc 
blackbody with a low characteristic temperature (~ 0.1 — 0.4 keV), 



plus a power-law continuum with a spectral index ~ 1.5 — 2.5 (e.g. 
|Feng & Kaaret 2005). The high luminosity and low temperature of 
the disc component were interpreted as supporting evidence for the 
existence of IMBHs (Miller et al.|2003||2004"l >. However, in the last 
few years, Chandra and XMM-Newton observations revealed new 
behaviours, showing the existence of X-ray spectra with rather pe- 
culiar properties. For some of the brightest ULXs, equally accept- 
able fits of their spectra may be obtained with (physically) differ- 
ent models, that suggest the presence of an optically thick corona, 
a fast ionized outflow o r a slim disc (e.g. I Stobbart et al.||2006| 
|Goncalves & Soria|20"06||Mizuno et al.|2007 1. Recently, |Gladstone| 
|et al.H2009) have shown that disc plus Comptonized corona models 
fit the data of the highest quality XMM-Newton spectra of several 
ULXs well, suggesting that the corona is cool and optically thick. 

In order to increase our understanding of ULXs and shed light 
on the mechanism at the origin of their powerful emission, it is cru- 
cial to investigate the evolution of their accretion flow through the 
variability of their X-ray spectra. In this work we present a sys- 
tematic analysis of the X-ray spectra of two ULXs in the galaxy 
NGC 1313, using three years of XMM-Newton observations. The 
barred spiral galaxy NGC 1313 hosts three ULXs, one of which 
(NGC 1313 X-3) is a known supernova (SN 1978K) interacting 
with the circumstellar medium. The other two ULXs, NGC 1313 
X-l and X-2 (X-l and X-2 hereafter), are located in different po- 
sitions, close to the nucleus (~ 50") X-l and at the outskirts of 
the host galaxy (~ 6') X-2. NGC 1313 has been observed several 
times over the years by XMM-Newton and a sufficient number of 
X-ray spectra are now available to attempt a characterization of the 
spectral variability of the sources hosted in it. 

|Feng & Kaaret| (2006] fitted a sequence of 12 XMM-Newton 



2 Fabio Pintore & Luca Zampieri 

Table 1. Log of the observations. 



No. 


Obs ID 


Date 


Exp" 
(ks) 


Instr. (X-l) b 


Instr. (X-2) 6 


Net counts(X-l) 


Net counts(X-2) 


1 


0150280101 


11/25/2003 


1 


M1/M2 


Ml/M2/pn 


651,637 


578, 643, 734 


2 


0150280301 


12/21/2003 


7.4 


pn 


Ml/M2/pn 


7648 


2960, 2874, 6304 


3 


0150280401 


12/23/2003 


3.2 


pn 


Ml/M2/pn 


2223 


2215,2151,3030 


4 


0150280501 


12/25/2003 


1.7 


pn 


Ml/M2/pn 


484 


1379, 1378,904 


5 


0150280601 


01/08/2004 


6.5 


pn 


Ml/M2/pn 


5524 


2714, 1694, 1690 


6 


0150281101 


01/16/2004 


2.7 


Ml/M2/pn 


Ml/M2/pn 


1770, 1993, 2366 


887,713, 1013 


7 


0205230201 


05/01/2004 


7.8 C 


M1/M2 


M1/M2 


1310, 1489 


700, 783 


8 


0205230301 


06/05/2004 


8.7 


Ml/M2/pn 


Ml/M2/pn 


5446,5151,8685 


3800, 3798, 819 


9 


0205230401 


08/23/2004 


3.8 


Ml/M2/pn 


Ml/M2/pn 


2397, 2643, 2323 


1252, 1316, 1173 


10 


0205230501 


1 1/23/2004 


12.5 


M1/M2 


Ml/M2/pn 


3013, 3179 


1523, 1594,4164 


11 


0205230601 


02/07/2005 


9.0 


Ml/M2/pn 


Ml/M2/pn 


5025,2435,2187 


3702,3697, 8224 


12 


0301860101 


03/06/2006 


17.2 


pn 


Ml/M2/pn 


3980 


4043,4374, 11426 


13 


0405090101 


10/16/2006 


78.6 


Ml/M2/pn 


pn 


24277,25141,58116 


53271 



a GTIofEPIC-pn 
pn = EPIC-pn camera; M1/M2 = EPIC-MOS1/MOS2 camera 
c GTI of EPIC-MOS 



observations of X-l and X-2 with a power-law plus multicolor disc 
blackbody model and found an anti-correlation between the lumi- 
nosity and the inner temperature of the MCD component. For this 
reason they concluded that such component does not originate in 
a standard accretion disc. The optical and X-ray variability of X-2 
was also investigated by Mucciarelli et al. (2007) to constrain the 
properties of the donor star and the binary system. They found that 
the power-law component hardens as the flux increases, opposite to 
what usually shown by Galactic BH XRBs. 

Another crucial issue related to ULX formation are the prop- 
erties of the environment in which ULXs are embedded. Claims re- 
garding the correlation of ULXs with low metallicity environments 
have been recently reported. Sw artz et al.| p008 1 found that within 
the Local Volume the specific ULX frequency decreases with host 
galaxy mass above ~ 1O 8 ' 5 M0 , meaning that smaller, lower metal- 
licity systems have more ULXs per unit mass than larger galaxies. 
|Mapelli et al] ( |2009^ and |Zampieri & Roberts]j2009^ suggested that 
at least a fraction of ULXs may be powered by massive stellar BHs 
formed from the direct collapse of low-metallicity massive stars. 
Using binary synthesis calculations, Linden et al. (2010) proposed 
another interpretation, in which the number, the lifetime and (to 
a less extent) the luminosity of high mass XRBs are enhanced at 
low metallicities. Only few measurements of the metallicity in the 
ULX environment are available, and results are not conclusive (see 
e.g. |Zampieri & R oberts 2009 and references therein). While opti- 
cal observations provide probably the best means to perform such 
measurements, also the abundances inferred from the detection of 
K-shell photoionization edges of intermediate mass or heavy ele- 
ments in the X-ray absorption spectrum of ULXs can be a viable 
tool, if high signal-to-noise spectra are available. This was already 
attempted by [Winter et al.| ( [2007| > for a sample of 14 ULXs with 
XMM-Newton spectra, obtaining values that match the solar abun- 
dance, but no further investigation, especially using the high reso- 
lution RGS spectra, has been attempted since then. 

In this work we try to characterize the spectral variability of 
X-l and X-2 using all the available XMM-Newton data. Recently, 
detailed spectral variability analyses were performed by Feng & 
|Kaaret| ( |2009) and |Kajava & Poutanen| (|2009|> on XMM -Newton 
and Chandra data of some ULXs, by |Dewangan et al.] ( |2010| > on 
the last two XMM-Newton observations of NGC 1 3 1 3 X- 1 , and by 



|Vierdayanti et al.| ( [20 1 0[ > on XMM-Newton and Swift data of Ho 
IX X-l. We will try to constrain also the metallicity of the absorb- 
ing gas towards X-l and X-2 using the RGS and EPIC spectra of 
the longest observation, and stacking together all the RGS observa- 
tions. Some preliminary results of this investigation were reported 
in |Pintore & Zam pieri (2011). Because of some improvements in 
the present analysis, the results and interpretation reported here su- 
persede our previous ones. The plan of the paper is the following. 
In §[2] we summarize the data selection and reduction procedures, 
while in §[5]and §[4]we present our results and discuss them in §|5] 



2 DATA REDUCTION 

We re-analyzed all the available XMM-Newton spectra of the two 
ULXs hosted in the spiral galaxy NGC 13 13, X-l and X-2, with ho- 
mogeneous criteria. The 17 observations span a time interval of six 
years, from 17 October 2000 to 16 October 2006, but three obser- 
vations were excluded because of high flares contamination. Also 
the observation of October 2000 was excluded from the analysis 
because the calibration before December 2000 may be incomplete. 
Data were reduced using SAS v. 9.0.0. EPIC-MOS and EPIC-pn 
spectra were extracted selecting the good time intervals with a 
background count rate not higher than 0.45 count s _1 in the en- 
ergy range 10 — 12 keV. We set 'FLAG=0' in order to exclude the 
events at the CCD edge and the bad pixels. Spectra were extracted 
from events with PATTERN sC 4 for EPIC-pn (which allows 
for single and double pixel events) and PATTERN ^ 12 for 
EPIC-MOS (which allows for single, double, triple and quadruple 
pixel events). We used 35" and 30" circular extraction regions for 
X-l and X-2, respectively. For the background we chose a 65" ex- 
traction circular region on the same CCD chip where the source is 
located. If the sources were on or near a CCD gap, no spectra were 
extracted. For X-2 the off-axis angle is quite similar for most of 
the observations (not much larger than 0.4' for all but three obser- 
vations). For X-l, 10 observations have off-axis angles that differ 
no more than 1', while the remaining three show a more significant 
variation. However, their spectra do not seem to present any pecu- 
liarity possibly associated to variations in the fraction of encircled 
energy. RGS spectra were extracted using the rgsproc task with the 
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Figure 1. Optical depth r versus temperature of the corona kT cor for the X-ray spectral fits of X-l (left) and X-2 (right) with a diskbb+comptt model. 
Observations that have a very thick corona are represented with black triangles while observations with a thick corona are represented as red squares. 



option spectrumbinning=lambda. In this way it turns out to be pos- 
sible to combine the spectra of different observations. The EPIC 
spectra were grouped with a minimum of 25 counts per channel, 
while RGS spectra with 50 counts per channel. 

All the spectral fits were performed using XSPEC v. 12.5.1. 
To improve the counting statistics, whenever possible, we fitted the 
EPIC-pn and EPIC-MOS spectra simultaneously. EPIC spectral fits 
were performed in the 0.3-10.0 keV energy range, while RGS fits 
were limited to the 0.45-2.0 keV energy range. For each instru- 
ment, a multiplicative constant was introduced to account for pos- 
sible residual differences in the instrument calibrations. The con- 
stant of the EPIC-pn data set was fixed equal to 1 , while the other 
two are allowed to vary. In general, the differences among the three 
instruments are not higher than 10%. 



3 X-RAY SPECTRAL FITS 

3.1 Comptonization plus multicolor blackbody disc 

We analyzed 13 out of the 17 XMM-Newton observations of X- 
1 and X-2 adopting an absorbed multicolor blackbody disc plus 
Comptonization component, modelled with diskbb+comptt ( |Mit-| 
|suda et al.|[l984l |Titarchuk|[l994l > in XSPEC. This model was 
successfully adopted in previous investigations to describe phe- 
nomenologically the spectra of ULXs ( Sto bbart et al. 12006] |Glad-| 
stone et al. 2009 , Feng & Kaaret 2009 ). Although several of our ob- 
servations show acceptable fits using only the comptt component, 
in some cases adding the diskbb component leads to a significant 
improvement in the fit. Therefore, in order to perform a comparison 
within the framework of a unique spectral model, here we assume 
the diskbb+comptt as reference model for all the observations and 
check its consistency on the basis of the variability patterns of its 
spectral parameters. The interstellar absorption was modelled with 
the tbabs model in XSPEC. We fixed the Galaxy column den sity 



along the line of sight at 3.9 x 10 cm l Dickey et al. 



1990) and 



added a free absorption component to model the local absorption 
near the source. 

The results of the spectral fits with the diskbb+comptt model 
for the 13 XMM-Newton observations of X-l and X-2 are reported 
in Table[2] We tied the temperature of the disc (Tdisc) to that of the 
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Energy (keV) 




1 2 

Energy (keV) 

Figure 2. X-ray spectra (energy X photons) of X-l (top) and X-2 (bottom) 
for observations that have a very thick corona (black) and a thick corona 
(red). In both cases the comparison is between observations #9 and #12, 
respectively. 
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Table 2. Best fitting spectral parameters of NGC1313 X-l and X-2 in different observations obtained with the absorbed diskbb+comptt model. 



NGC1313X-1 


No. 


Date 


N H a 
(10 21 cm 2 ) 


(keV) 


kT d 

<v ± COT 

(keV) 


r e L x [0.3-10 keV] f L diBC [0.3-10 keV]f 
(lO^ergs- 1 ) (10 39 ergs- 1 ) 


X 2 /dof 



1 


1 1/25/2003 


2+0.2 

U - Z -0.2 


o.8oi° : °i 


9.03+ 400 


16.2+ 2 


1Q +0.3 


c nl.2 

a.u 4 . 


21.33/38 


2 


12/21/2003 


1 9+ 01 


0-199™ 


1 OO+0.02 

J - oo -0.02 


6 7+ 01 
"■'-0.1 


11 5+ - 9 
J-l-O-0.9 


n+3e36 
U 


195.47/226 


3 


12/23/2003 


4 0+ 02 


929+ - 004 
u.zza_ 004 


2 4 i+0.08 
z '* A -0.08 


6 2+ 03 
°- z -0.3 


14.2+22 


4 6 ' 8 
4 - D 0.9- 


64.38/75 


4 


12/25/2003 


4 3+0-4 


Q 179 +0.005 
Uil ' -0.006 


9 7+O.2 
Z ' ' -0.2 


6 2+ - 7 
"•-0.7 


14 o+5.'9 


6.6 2 ' 4 


12.88/12 


5 


01/08/2004 


n 7+0.1 
z - ' -0.1 


244+ - 003 
u ' z * -0.003 


45 +o.i 


4.3i° j 




9 7+0.4 
Z ''-0.3 


184.56/176 


6 


01/16/2004 


1 3+0.1 
0.1 


226+ 006 


9 9 o+0. 04 
Z - Z »_0.Q4 




9-3l : 8 6 


nl.6e38 
U 


194.71/197 


7 


05/01/2004 


q 9+O.2 
°' z -0.2 


216+ 003 
U ' Z1D -0.003 


4g +0.2 


4 8+ - 3 


8.0+i;« 


z -^0.4 


100.97/91 


8 


06/05/2004 


i o 7 +0.07 


242+0-003 
u,z *-0.003 


1 70+O.O2 


7 1+0.1 

'•• L -0.1 


15.0+1 


1 3 ' 2 


511.38/543 


9 


08/23/2004 


9 q+0.1 

z -°-o.i 


u.roo_ 002 


2 =Q+0.03 
z - d3 -0.03 


3 7 +0.07 
°- ,u -0.07 


5 4+ 06 


nl.4e38 
U 


233.91/232 


10 


1 1/23/2004 


9 9+0.1 

z - z -o.i 


3 67 +o.oo3 
u -°°' -0.003 


1 48+ - 04 
1 -* o -0.04 


14 6 +0 ' 8 


8 - 2 -0.9 


3 7O.3 


165.41/184 


11 


02/07/2005 


9 O1+0-09 

z - a -o.os 


n 9ns+o oo 2 
u.2U»_ 002 


9 71 +0.05 
z -' 1 -0.05 


7 0+ - 2 

'■ u -0.2 


8 9+ 1 


9 7+O.3 
z -'-0.3 


313.10/305 


12 


03/06/2006 


9 q+0.1 
z,o -0.1 


45O+0-004 

u.ioy_ 004 


1 01+ 003 
1 - ul -0.03 


15+* 7 


4 6+°' 6 


z -°-0.2 


168.29/139 


13 


10/16/2006 


n fi7 +0.002 
z -°'-0.002 


n 999+O.OOI 
U - ZZZ -0.001 


2 14 +0.0l 
z - 14 -o.oi 


8 44+0-06 
^ -0.06 


7 1+0.2 
'• x -0.2 


2 1+ 01 

z - i -o.i 


1490.95/1420 













NGC1313X-2 








No. 


Date 


N H a 




kT c 
cor 




L x [0.3-10 keV] e 


L dlsc [0.3-10 keV]/ 








(10 21 cm 2 ) 


(keV) 


(keV) 




(10 39 ergs- 1 ) 


(10 38 ergs" 1 ) 




1 


1 1/25/2003 


6+ 02 
u.o_ Q 2 


44+0-01 
u -^-0.01 


, qC+0.05 


10.8+° 7 


5 4+1.4 

M -1.2 


12l| 


58.72/62 


2 


12/21/2003 


1 6+ 01 


n qcn+O- 004 
U.o8U_ Q Q04 


1 5fi+°- 02 
i - OD -0.02 


ILOtp 


7 q+0.7 
'•°-0.6 


14+2 


386.05/395 


3 


12/23/2003 


5+ - 1 


27+ - 01 

u - z '-o.oi 


1 74+0.03 


io-olS:2 


77+1 
'•'-0.8 


n 6.4e+37 
U 


245.35/248 


4 


12/25/2003 


1 9+0.2 
^-0.2 


469+0 06 


3.02+°;* 


6 3+ 06 


q 7+O.8 
°-'-0.6 


14.4+ 2 


99.68/120 


5 


01/08/2004 


1 4+ - 1 


266+ 005 
U - ZDD -0.005 


1 89+ - 04 

1 - oy -0.04 


6 6+ 02 

°'°-0.2 


9 H+0- 4 
Z - S -0.4 


4 5+ 1 


175.27/201 


6 


01/16/2004 


9 1+0.2 
z - 1 -0.2 


0-186±g;S8g 


5.93+0' 2 


q q+0.2 

o.o_ Q 2 


q c+0.9 


4 7+2-3 
*-'-1.7 


84.94/85 


7 


05/01/2004 


1 0+ 03 
X - U -0.3 


0.381°$ 


1 44+0-07 
— 0.07 


7 9+0.7 
' - Z -0.7 


1 o+0.5 
l- S -0.4 


6 3+ 17 


55.41/46 


8 


06/05/2004 


1 1+ 01 


n 9q«+0.006 
U.238_ Q 00{j 


1 48+ ' 01 
J -*°— 0.01 


11 0+ ' 2 


n +0.7 


n +4.8e39 
U 


514.96/507 


9 


08/23/2004 


3 Q+0.1 


n 9n o+0.002 
U.2U»_ Q 0()2 


1 96+ 05 


7 4 +0.3 


9 S+0- 4 
z - 8 -0.4 


121 


150.72/131 


10 


1 1/23/2004 


9 q+0.1 
z -°-0.09 


226+ ' 002 


3 84+ - 08 
°-° -0.08 


44+0-1 


2 6 0+ - 3 
z -°-0.3 




210.89/237 


11 


02/07/2005 


1 69+ - 09 

±,u -0.08 


goo+0.004 
u.uoo_ 0Q4 


, =,-+0.02 
i - OO -0.02 


13 0+ ' 5 


a i+0.4 


28+ 2 


490.92/486 


12 


03/06/2006 


1 54+0-08 


n 9CQ+0.006 
U.2o»_ 006 


1 42+0- 01 
1 -^ z -0.01 


11 5+0-2 
41 -°-0.2 


7 c+0.6 
'■°-0.6 


5 1+ 1 


577.03/603 


13 


10/16/2006 


1 89+ - 04 
1 - oy -0.04 


Ogi+0.02 
u - D1 -0.O2 


1 61+ 001 


19 q+0.3 
lz -°-0.3 


7 9+O.2 
' - Z -0.3 


28lo:2 


841.83/854 



a Column density; b Inner disc temperature; c The seed photons temperature Tq is assumed to be equal to T c n sc . d Temperature of the corona; e Optical 
depth of the corona; f Unabsorbed total x-ray luminosity in the 0.3 -10 keV range; 9 Unabsorbed disc luminosity in the 0.3-10 keV range. 
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Figure 3. Optical depth r versus unabsorbed total luminosity (in the 0.3-10 
keV range) for X-2. The black points represent the very thick corona state 
while the red crosses are the thick corona state. 



seed photons (To) for comptonization. We will comment on this 
choice at the end of this Section. The reported errors are at the 90% 
confidence level for one interesting parameter. We noted that there 
are often several local minima with close values of the x 2 > some- 
times with evidence for both a strong/warm and a weak/cool (or no) 
disc. After a careful inspection of the \ 2 surface, we found the ab- 
solute minima for each observation reported in Table [2] However, 
one should be aware that the actual uncertainty on the disc parame- 
ters caused by the topology of the % 2 surface is at least a few times 
larger than the formal error reported in the Table. We tried to adopt 
also more physical models, such as the eqtherm (Coppi 2001 1 and 
the DKBBFTH Pone & Kubota|2006| >. We found that the count- 
ing statistics of most of the spectra is inadequate to constrain the 
parameters of eqtherm and DKBBFTH (see also below). 

The longest observation of X-l has a short segment of ~ 3.5 
ks within the first 13 ks that shows some intrinsic variability (a 
slight systematic hardening at high energies) and was then removed 
from the analysis. We note also that observation #1 has a very low 
value of Nh compared to the other observations and a low counting 
statistics. It converges towards a local minimum with a very strong 
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Figure 4. X-ray unabsorbed luminosity evaluated in the 0.3-10 keV energy band (a distance of 3.7 Mpc was assumed; Tully 1988 1 for all the XMM-Newton 
observations of X-l (up-left) and X-2 (up-right). Observations with a very thick corona are represented with black triangles while observations with thick 
corona are represented as red squares. All the (unabsorbed) luminosities are evaluated in the 0.3-10 keV energy band. 



disc component and with parameters unlike those of all the other 
spectra. For this reason, we excluded it from the following analysis. 

Figure [T] shows a plot of the optical depth r versus the tem- 
perature of the corona kT cor obtained from the best fits of both 
sources. There appears to be well defined locii in the kT cor -T plane, 
indicating the existence of a somewhat ordered behaviour in the 
spectral variability of the corona. The spectra of X-2 populate two 
distinct regions characterized by very large optical depths (r > 10) 
and low temperatures (kT cor ~ 1.5 keV), on one side, and smaller 
optical depths (r < 8) and a range of temperatures (kT cor ~ 1.5- 
6 keV), on the other side. In both cases the corona turns out to 
be optically thick. In the following we refer to these two regions 
as "very-thick" and "thick" corona states, respectively. X-l may 
also be interpreted within the same framework. An analysis of the 
X-ray flux variability shows that both sources tend to be signifi- 
cantly more variable when the coronae are thicker (fractional vari- 
ability >9-10%, see table[3j, which would strengthen the analogies 
between their behaviours. However, the clustering of the observa- 
tions in the kT cor -r plane of X-l is not so well defined and the 
region with high optical depths is populated by only two observa- 
tions. Furthermore, although not included in Figure [T] because of 
the low counting statistics, our interpretation is not consistent with 
the spectral paremeters inferred from observation #1. So, present 
data are not sufficient to establish whether X-l has indeed a bi- 
modal behaviour similar to that shown by X-2 or the changes in 
the optical depth should be interpreted in a different way (e.g. as 
due to different variability patterns or accretion geometries; |Feng| 
|& Kaaret||2006| |Dewangan et al.||2010| >. In the following we will 
continue to distinguish between the observations of X-l with a very 
thick corona (r > 10) and a thick corona (r < 8), being aware that 
they may represent physically different states with respect to those 
of X-2. 

An example of the spectral shapes when the two sources are in 
different positions on the kT cor -r plane is shown in Figure|2] Again 
there are analogies and differences between X-l and X-2. In both 
cases the spectra for very high coronal depths are bell-shaped, with 
a clear turn-over at > 3 — 4 keV (e.g. Stobbart et al. 2006), whereas 
the spectra in the low-r region are steeper and do not show strong 
evidence of curvature at high energies. Only for X-2 the two spec- 
tral states appear to correlate with total luminosity, the very-thick 



corona state being more luminous (see Figure[3]l. In X-l there is no 
significant dependence of the spectral shape on Lx- This is evident 
also from the behaviour of the spectra shown in Figure|2] While for 
X-2 the thick corona spectra stay always below those in the very 
thick corona state, for X-l there is a sort of crossing/pivoting point 
of the observed spectra (see also Kajava & Poutanen 2009) at ~ 2 
keV. Therefore, the total counts in the thick corona state of X-2 
are clearly smaller than those in the very thick corona state, while 
in X-l the deficit of photons observed at low energies in spectra 
with larger optical depths is compensated by the excess of photons 
at high energies. Figure [4] shows the light curves of X-l and X-2 
computed from the best fitting diskbb+comptt model (a distance 
of 3.7 Mpc was assumed; [Tu"lly|1988| >. On average, X-l has a lu- 
minosity (~ 10 40 erg s" 1 ) higher than that of X-2 (~5x 10 39 
erg s~ ). Variability of a factor ~ 3 and ~ 5 is observed for X-l 
and X-2, respectively. We further investigated the behaviour of the 
soft component in both sources. The discs are soft or warm, with 
temperatures of ~ 0.2 - 0.5 keV for X-l and ~ 0.2 - 0.6 keV 
for X-2. Three spectra of X-l are consistent with zero normaliza- 
tion (or absence) of the soft component. In eight observations the 
luminosity of the soft component is ~ 10 39 erg s _1 and it repre- 
sents a significant fraction of the total flux (> 30%). Excluding the 
spectral fits that return zero normalization of the soft component, 
there is no evidence of correlation or anti-correlation between the 
disc (or total) luminosity and the inner disc temperature. In X-2 the 
luminosity of the soft component is a significant fraction (> 30%) 
of the total luminosity in six observations. Two observations have 
normalization consistent with zero. Excluding them, the disc lumi- 
nosity appears to show a weak power-law correlation with the inner 
temperature, Ldisc oc T^:^ ' 3 , However, the correlation is uncer- 
tain, as the two observations that substantiate it (those with higher 
kTdisc) have also rather shallow minima in x 2 > admitting both a 
strong/warm and a weak/cool disc fit with close values of the \ 2 - 
No correlation is found using the total luminosity. 

3.2 Effects of varying the ratio of seed photons temperature 
to the disc temperature 

The comptonizing coronae of the comptt model turn out to be op- 
tically thick. This poses a problem, as in these physical conditions 
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Table 3. Fractional variability of NGC 1 3 1 3 X- 1 and X-2. 



Very Thick corona 


Energy band" 


X-2 (Obs. 13) 
counts 


F var (per cent) 6 


X-l (Obs. 12) c 
counts s~ 1 


F V ar (per cent) 6 


0.2-10 keV 


0.978 ± 0.005 


9.77 ±0.07 


1.47 ± 0.09 


14.6 ±0.6 


Thick corona 


Energy band" 


X-2 (Obs. 10) 
counts 


-F„ a r(percent) 6 


X-l (Obs. 13) c ' e 
counts s~ 1 


F V ar (per cent) 6 


0.2-10 keV 


0.377 ± 0.006 


< 9.3 d 


0.77 ± 0.07 


9 O+0.9 
z -°-1.9 



a The adopted energy band is for a direc t comp arison with |Dewangan et al.|j2010) 6 Calculated from the background subtracted EPIC-pn light curve, with 
200 s time bins. c From |Dewangan etal.^OlO} . d 3 o upper limit. e Fractional variability computed for the whole observation, including the short segment 

of 3.5 ks that we left out from the spectral analysis. 



the disc underneath the corona is masked by it. Therefore, the tem- 
perature of the disc component refers to the outer visible part of the 
disc, while the temperature of the seed photons To is not necessar- 
ily equal to Tdi sc - For this reason, Gladstone et al. (2009) adopted 
the DKBBFTH model in which the corona is assumed to cover the 
inner disc. We tried to apply the same model to our spectral se- 
quence, but found that it converges to plausible physical values of 
the parameters only for the highest quality spectra. We then tried to 
repeat our analysis with the diskbb+comptt model and disconnect- 
ing the two temperatures, but this leads to difficulties in finding a 
global minimum (see also the analysis of NGC 5204 X-l in Feng & 
Kaar et|2009| >. Finally, we decided to test our results tying the two 
temperatures with a fixed proportionality constant. 

If the corona is optically thick and is absorbing a constant 
fraction / of the accretion power, the actual inner disc tempera- 
ture Tj is lower than the inner temperature of th e disc in absence 
of the corona Ti, as t[ = Ti(l - /) 1/4 (e.g. [Gladstone et al. 
2009). If Ri and Rt are the inner disc radius and the trunca- 
tion radius of the corona respectively, from the relation T(R) « 
Ti (R/ Ri)~ 3 ' 4 , valid for a disc in absence of the corona, we obtain 

Ti » T t (Rt/Ri) 3/4 - Thus, we have T x « T t (Rt/Ri) 3/4 (1 - 
J) 1 / 4 . Assuming that / < 90% and that the corona is compact 
(Rt < 4i?i), the inner disc temperature Ti , which is also the seed 
photons temperature, is larger than the temperature of the outer vis- 
ible disc It, and 1 < T x /Tt = T /T dt3C < 2.5. We then at- 
tempted to perform some additional fits with two fixed values (1.5, 
2) of the ratio 8 = To /Tdisc as representative of this situation. 

We found that the spectral fits with a diskbb+comptt model 
and 8 = 1.5, 2 are statistically acceptable and the coronae are still 
optically thick, with the inferred optical depth weakly depending 
on To /Tdisc- The parameters of the corona change by no more 
than ~ 60%, with typical variations of < 20 — 30%. However, 
the dependence of spectral states on the total luminosity that char- 
acterizes the behaviour of X-2 for equal temperatures is lost. For 
X-l, the temperature of the disc continues to not correlate with the 
luminosity, as found for equal temperatures. For X-2, varying 8 the 
very-thick state becomes less populated, as some observations pre- 
viuosly in that state move to the thick corona state. For 8 — 1.5 
there is a weak correlation between disc luminosity and tempera- 
ture (L d isc oc T^ 4 * 0,4 ), while if 8 = 2 the correlation disappears. 
At variance with the correlation found for equal temperatures, the 
correlation for 8 — 1.5 is not critically dependent only on two ob- 
servations. 




Energy (keV) 

Figure 5. Stacked RGS 1 spectrum of X- 1 , along with its best fitting model 
(see text for details). 



These results do not change significantly for slightly more ex- 
tended coronae. On the other hand, if the corona is very extended 
and optically thick, the value of 8 may become so large that the 
disc component falls essentially outside the XMM-Newton band- 
pass (unless the disc-corona coupling is very strong). In these con- 
ditions some spectra are no longer well fitted by the model. 



4 CHEMICAL ABUNDANCE ESTIMATES 

We analyzed the EPIC-pn and RGS data of all the XMM-Newton 
observations of X-l and X-2 in an attempt to use them for deter- 
mining the chemical abundances in the local source environment. 
Following [Winter et al.|p007] >, for the EPIC-pn 5000 and 40000 
counts are necessary to observe the Oxygen K-shell and Iron L- 
shell absorption edges at 0.538 keV and 0.851 keV, respectively. 
Observation #13 has ~ 58000 (53000) EPIC-pn counts for X-l 
(X-2), and hence we can perform a chemical abundance analysis 
on the EPIC-pn spectrum similar to that presented in Wint er et al.| 
( |2007| l. In the spectral fits the tbabs absorption model is replaced 
with tbvarabs that allows to vary the chemical abundances (and 
grain composition). We set alternatively the abundance of Oxygen 
or Iron to zero. The spectrum was then fitted with the EPIC contin- 
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Table 4. Abundances inferred from the Oxygen K-shell photoionization 
edge (0.538 keV) from the X-ray spectral fits of NGC 1313 X-l and X-2. 



NGC 1313 X-l 






12 + \og(0/H) b 


RGS (observation 13) 
EPIC (observation 13) 
RGS (stacked) 


U -'°-0.25 
u -'°-0.03 
U.DO_ 23 


S 7+0-1 

o vo +0.02 
»-'^_0.02 

8 -'-0.2 


NGC 1313 X-2 






12 + log(0/H) 


EPIC (observation 13) 


46+ - 05 


o fi/ ,+0.06 
S - b4 -0.05 



a Absorption depth at the threshold energy. 

6 For the solar abundance we assume 12 + \og(0/H) = 8.92 (Zq = 
0.02). 



uum model in Table [2] (keeping all parameters, but normalizations, 
fixed) plus an absorption edge, that accounts for the observed ab- 
sorption feature. The parameters of the edge are then used to com- 
pute the abundance. 

Using this approach, we found that the O and Fe abundances 
inferred from the EPIC-pn spectrum of X-l are consistent with a 
sub-solar metallicity environment (Z ~ 0.6Zq from the O edge; 
Table EJ. The \ 2 of the fit with the absorption edge is 1488 for 
1423 dof, while that without the edge is 1759 for 1424 dof. For 
X-2 we find a subsolar abundance for Oxygen (Z ~ 0.5Z@; see 
again TableQ and an Iron abundance consistent with zero. The fits 
returns x 2 = 839 (857 dof) with the O edge and X 2 = 953 (858 
dof) without it. 

We tried to analyze also the RGS data of X-l and X-2 using 
the same technique. As the sources are faint, the RGS net count 
rate is quite low (~ 1.5 x 10~ 2 count s" 1 ). Therefore, only the 
last observation (observation #13, 122 ks) reached a reasonable 
counting statistics in the RGS for the brightest ULX (X-l), while 
no useful analysis could be performed on X-2. We tentatively iden- 
tified two features in absorption in the RGS spectrum of X- 1 , as- 
sociated to O I (0.535 keV) and Fe I (0.709 keV). There may be 
also two other lines in emission at the characteristic energy of O 
VIII K a (0.653 keV) and Si K a (1.748 keV), but their significance 
is very low. From the Oxygen edge we found an abundance below 
solar (Z ~ O.GZq) for the absorbing material towards X-l (see Ta- 
ble |4j, although the statistical improvement obtained including the 
absorption edge is small. 

We tried to improve the analysis stacking together the RGS 
spectra in such a way to increase the counting statistics. We used 
only the RGS1 spectra, because our best diagnostic is the Oxy- 
gen K-shell photoionization edge at 0.538 keV and the response 
of the RGS2 has some problems precisely in this energy range 
(0.5-0.6 keV). The stacked spectra of X-2 have not enough count- 
ing statistics for a meaningful analysis. For X-l, we combined 
spectra together using the rgscombine command, that appropri- 
ately accounts for the response matrices and backgrounds of the 
different observations. We had to exclude several observations 
(#2, 3, 4, 5, 6, 7, 11, 12) that caused technical problems in the pro- 
cessing with XSPEC. Observation #1 was also excluded for the 
reason explained in the previous Section. The final combined RGS 
spectrum of X-l has a total exposure time of 168 ks and ~ 2861 net 
counts. We fitted it with a diskbb+comptt model with parameters 
fixed and equal to the mean values obtained from the EPIC spectral 



fits of X-l, setting the Oxygen abundance to zero and adding an 
absorption edge (Figure 13}. We found that the Oxygen abundance 
is again subsolar (Tablepk. The \ 2 of the fit with the absorption 
edge is 139 for 103 dof, while that without the edge is 146 for 104 
dof. However, the actual abundance is rather sensitive to the tem- 
perature of the soft component and its uncertainty becomes large if 
kTdisc is included in the spectral fit. 



5 DISCUSSION 

The analysis of all the XMM-Newton observations shows that the 
spectra of X-l and X-2 can be well reproduced by a Comptoniza- 
tion model plus a soft disc component in which the coronae are 
always optically thick. Both sources appear to show a well de- 
fined behaviour in the optical depth (r) versus corona temperature 
(kT cor ) plane. For X-2 we clearly identified two states that charac- 
terize the spectral variability and appear to have also different short 
term variability properties: a "very-thick" (more variable) corona 
state in which r > 10 and kT cor ~ 1.5 keV, and a "thick-corona" 
state in which r < 8 and kT cor ~ 1.5 — 6 keV. We note that 
a morphological classification in terms of a harder brighter state 
and a softer dimmer state, based on a single power-law fit, has al- 
ready been proposed for X-2 (Feng & Kaaret 2006 1. Here we offer 
a physical explanation of these two states in terms of varying pa- 
rameters of an optically thick corona. The behaviour of X- 1 may be 
interpreted within the same framework but, with presently available 
data, the observed changes in the coronal optical depth could be ex- 
plained also in a different way (e.g. as due to different variability 
patterns or accretion geometries; Feng & Kaaret 2006; Dewangan 
|et al.|2010| >. 

For X-2, which is on average less luminous than X-l, the two 
spectral states of the corona appear to correlate with luminosity. 
The optical depth of the corona increases as the total luminosity 
goes up, as expected if the corona responds rapidly to an increment 
in the instantaneous accretion rate. The behaviour of the disc com- 
ponent is also different in the two sources. While in X-l the lumi- 
nosity and temperature of this component do not correlate, for X-2 
we find that Ldi SC <x T^r^r ' 3 . However, as already mentioned, 
this result is uncertain as the two spectra that substantiate the cor- 
relation admit both a strong/warm and a weak/cool disc fit with 
close values of the % 2 ■ 

A spectral analysis based on the assumption that the seed pho- 
ton temperature (To) is equal to the disc temperature (Tdi 3C ) has 
some inconsistencies because, if the corona is optically thick, the 
innermost part of the accretion disc is actually not visible and, 
hence, it is not necessarily true that To = Tdisc- For this reason, 
we repeated our analysis with the diskbb+comptt model tying the 
two temperatures with a proportionality constant, that was fixed 
assuming a compact corona energetically coupled to the disc (dis- 
connecting the two temperatures leads to difficulties in finding a 
global minimum because of the low counting statistics of several 
observations). We find changes in the corona parameters of < 60%, 
with typical variations of < 20 — 30%. However, the dependence 
of spectral states on the total luminosity that characterizes the be- 
haviour of X-2 for equal temperatures is lost. While for X-l the 
temperature of the disc continues not to correlate with the luminos- 
ity, for X-2 the correlation persists up to values of To /Tdi SC < 1.5, 
consistent with strong disc-corona coupling (/ ~ 90%). 

We emphasize that the model adopted in the present investi- 
gation is not entirely physically consistent. Besides the issues of 
the relation between To and Tdisc that we tried to address as de- 
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scribed above, there may be other caveats, such as the input pho- 
tons not being a mono-temperature Wien distribution (as assumed 
in comptt), or the disc structure being different from a standard one 
at high accretion rates, or the origin and precise location of the 
thick corona/wind being unkown. However, our purpose was not 
to adopt the most physically consistent spectral model for ULXs, 
which does not exists yet, but a simplified one that can reflect the 
underlying physics and can be tested on the basis of the observed 
spectral variability patterns, limited by the current data quality. 

A full explanation of the observed spectral variability patterns 
of X-l and X-2 is beyond the scope of the present investigation 
and appears not easy. Here we simply propose some interpretations 



within the framework of recent work on the subject (e.g. |Feng & 



Kaaret|2009l [Gladstone et al.|2009| [Dewangan et al. 2010| |Vier- 



dayanti et al |2010f . X-l shows higher average isotropic luminos 
ity (~ 10 4U erg s -i ) and smaller variability. The different spectral 
shapes do not correlate with luminosity, and the temperature and lu- 
minosity of the soft component do not vary together. These findings 
are consistent with X-l being in the ultraluminous regime JGlad-| 
|stone et al.|20 09l, accreting at super-Eddington rates and launching 
powerful winds from a disc embedded in an optically thick corona. 
The accretion disc may be there, but in a different physical regime 
and covered by a wind/corona, while the soft component may ac- 
tually represent emission from the wind itself. The larger X-ray 
fractional variability observed when the wind/corona is optically 
thicker may be consistent with increasing obscuration of the source 
caused by the turbolent wind itself at higher accretion rates ( |Mid-| 
|dleton et al.|20ri> ,b). 

As far as X-2 is concerned, it behaves in a way similar to X- 1 
on the T-kTc 0r plane and also the fractional X-ray variability of the 
two sources appears to be comparable. Hence the corona may be in 
a similar physical state. Again, this may be consistent with the pic- 
ture in which we are seeing the two sources under a high inclination 
angle and the temporal variability may be produced by blobs in the 
wind that intersect our line of sight to the hottest central regions of 
the source (Middleton et al. 201 1 i,b). However, as noted above, X- 
2 has a lower average luminosity and, for To < 1.5Tdisc, it shows 
a weak correlation between the luminosity and inner temperature 
of the disc component. This may be interpreted as X-2 being in 
a less extreme regime, in which the average accretion rate is at 
around or slightly above the Eddington limit. The soft spectral com- 
ponent may then physically represent a true accretion disc which 
is partly visible and has a characteristic temperature in the range 
0.2-0.6 keV. Clearly, some time-dependent interaction between the 
disc and the corona is likely to occur (e. g. a slight expansion on a 
dynamical timescale of the corona as the accretion rate increases) 
and may be responsible for the observed slope of the temperature- 
luminosity relation (Ldisc oc Tjif c ), which is different from that of 
a standard disc with fixed inner radius (Ldisc oc Td isc ). Clearly, as 
the corona is always optically thick, it is not possible to use the disc 
parameters for estimating the BH mass. 

We note that the disc temperature-luminosity correlation in X- 
2 disappears using a simple power-law to describe the spectrum of 
the corona, as found by |Feng & Kaar et (2006) (and by |Feng "&] 
|Kaaret| (2009"1 > for IC 342 X-l), probably because the power-law 
does not take into account the spectral curvature at high energies. 
The correlation of the optical depth to the corona with luminosity 
that we found for X-2 (assuming To = Tdi SC ) is in agreement with 
the results obtained for Ho IX X-l by Vierday anti et al.| ( |2010| >. 
Also IC 342 X-l shows some hint of an increase in the coronal 
depth as the X-ray luminosity increases < |Feng & Kaaret1|2009| >. 
These similarities in spectral variability patterns appear to be con- 



sistent with the proximity of these three ULXs in the spectral se- 
quence proposed by Glads tone et al.| (|2009 I. 

We used the RGS high spectral resolution to attempt an esti- 
mate of the metallicities of the local environments of X-l and X-2. 
Because of the low signal-to-noise ratio, only the last, longest ob- 
servation of X-l could be used. The analysis was performed also 
on the EPIC spectrum. The metallicity in the X-l environment is 
consistent with being below solar. The last EPIC spectrum of X- 
2 was also analyzed with the same method and suggests subsolar 
metallicity. In an attempt to increase the significance of our results 
we also stacked together the RGS spectra of some observations of 
X-l and performed measurements of the metallicity on the stacked 
spectrum. Also in this case the Oxygen abundance turns out to be 
below solar, although it is rather sensitive to the temperature of the 
soft component. Our estimates are in agreement with the abundance 
measurements from HII regions in NGC 1313 that give values that 
are all subsolar (Pilyugin 2001 ; Hadfield et al. 2007; Ripamonti et 
|al.|2010| >, but are smaller than the slightly supersolar metallicity of 
the X-l enviroment found by W inter et al.| ( |2007] >. The difference 
with the latter Authors may be due to the fact that we analyzed an 
observation with higher counting statistics and adopted a different 
spectral model. 

A couple of ULXs show periodic intensity variations in X-rays 
which are considered as signatures of the orbital period {Kaaret et] 
|al.|2006a|b| [Kaaret & Feng|2007| [Strohmayer|2009^ . X-ray spec- 
tra seem to change regularly with orbital phase. It is possible that 
these phase related variations may affect any result obtained from 
snapshot observations such as those presented here. Furthermore, 
as mentioned above, the observed spectral variability patterns of 
X-l and X-2 are rather complex and their interpretation does not 
appear directly comparable to that of Galactic XRBs. For these 
reasons, the acquisition of new, high quality spectra through a ded- 
icated X-ray monitoring programme is definitely needed. 
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